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ABSTRACT 

We demostrate that the measurement of fluctuations of polarization due to the galactic interstellar 
scintillations may be used to study the structure of the radiation field at compact radio sources. We 
develop a mathematical formalism and demonstrate it on a simple analytical model in which the scale of 
the polarization variation through the source is comparable to the source size. The predicted amplitude 
of modulation of the polarized radiation flux is ~ 20% x tTs x rUsc , where tTs is the characteristic degree 
of polarization of radiation at the source and ttIsc is the typical modulation index due to scattering, i.e., 
7TT.sc — 1 for diffractive scintillations and mgc < 1 for refractive scintillations. 
Subject headings: ISM: general — scattering — technique: polarimetric — technique: interferometric 



1. INTRODUCTION 

Interstellar scintillations of compact radio sources, 
which w ere first identified in the sh ort-term variability of 
pulsars ( Scheuer 1968| ; Ricket 1969| ) and explained as ran- 
dom scattering and diffraction from the electron density 
fiuctuations in the ionized interstellar plasma, have been 
used in studies of galac tic and extra-galactic objects (see. 
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e.g., Cordes et al. 1985; 1986 ), mostly pulsars, and also ac 



live galactic nuclei (see, e.g., phosh fc Gopal-Krishna 1990 
and references therein), some comp act extra- galactic radio 
sources (e.g., [Rickctt et al. 1995), 7-ray burst remnants 



( [Goodman 1997| [Waxman ct al. 1998), and th e accretion 
disk around the galactic center ( |Lo et al. 1998| ). However, 
this technique is (mostly) limited to the measurements of 
the intensity variations, while the information contained in 
the polarization properties is usually lost. In this letter we 
extend the interstellar scintillation observation technique 
to include the detection of polarization fluctuations, which 
will allow to study the structure of the radiation field at 
the source. For synchrotron radiation, for instance, this 
allows to infer the geometry and correlation properties of 
the magnetic fields, — the information which is difficult 
to obtain by other techniques becau se of resolution limits 
(see, e.g., Medvcdev & Loeb 1999 for an application to 
7-ray burstejl 

In this letter we develop a formalism of the interstellar 
scintillations of polarization. For simplicity, we do not con- 
sider frequency-related effects. Using a simple model of the 
source radiation, we analytically calculate modulation in- 
dices. We demonstrate that the amplitude of fluctuations 
of polarization is maximum when the angular correlation 
scale of radiation at the source (which is often the source 
size) is comparable to the typical scale of 'speckles' pro- 
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duced at the observing plane by the interstellar scattering, 
and amounts up to ~ 20% x tTs x rUsc modulation, where 
TTs is the typical degree of polarization at the source and 
rrisc is the typical modulation amplitude which is ~ 1 for 
diffractive scintillations and usually ~ 0.1 for refractive 
scintillations. 

The rest of this letter is organized as follows. In g2 we 
present a mathematical technique, which we then demon- 
strate on a simple analytical model in §0. We present our 
conclusions in §0. 

2. POLARIZATION SCINTILLATION FORMALISM 

The interstellar scintillations arise when fluctuations in 
the electron density randomly modulate the refractive in- 
dex of a medium. As a result of diffraction (diffractive 
scintillations) and random focusing (refractive scintilla- 
tions) of electromagnetic waves, a point source produces a 
pattern in an observing plane which consists of randomly 
located bright and dim spots — speckles — and observed 
as fluctuations of the source brightness when an observer 
moves through the pattern. The characteristic correlation 
length of the pattern, ro(A), is determined by the regime of 
scattering, statistical properties of the interstellar turbu- 
lence, and the obser ving wavelength , A (for more details, 
see a nice review by Narayan 1992 ).^ It is convenient to 
introduce the 'angular' size of a speckle, 9^ — r^/D, where 
D is the distance to the scattering screen, and we define 
^s to be the angular size of a source. 

Let us assume that the emission is polarized, and the 
sense of polarization varies through the source, so that the 
net (average) polarization vanishes.^ If 0s ^ ^07 the source 
is effectively point-like. In this regime the intensity vari- 
ations are large, while the polarization signal is weak due 
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^Let us introduce the Fresnel scale, rp = {\D/2tt) ' (where D is the distance to the scattering screen) and the diffractive length, r^ifi, 
which represents the transverse separation for which the root mean square phase difference (produced by scattering) is equal to one radian. In 
the weak scattering regime (r;jiff ^ rp), the spec kle size is of order of the Fre snel scale, ro ~ rp. In the strong scattering regime (r^ifj <g; rp), 

For the short-term diffractive scintillations, the speckle size is 



Joodman & Narayan 1985) 
It is ro 



ref = ^p/'"diff ^ ^diff- Note that the apparent angular size of the 



the diffractive and refractive effects are distinct ( ^ 

''0 ~ ''diffi while for the longer-term refractive scintillations, 

scattered point source image is 6scatt ~ r^^f/D for both diffractive and refractive scintillations. 

^If the net degree of polarization is nonzero, one can always subtract the constant contribution. The variation of this polarized contribution 
is identical to the total flux variation. 
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to averaging. In contrast, when 9s ^, Op, different elements 1995). In calculating the modulation indices below, we 



in the source produce identical patterns, which are ofF- 
set by D 59, according to the element angular separation, 
69. The overall pattern is, thus, obtained from the super- 
position (convolution) of individual, incoherent patterns. 
An observer will measure fluctuations of the direction of 
polarization, while the intensity exhibits less modulation. 
When 9s ^ 9q, the offsets become much larger than the 
speckle size, so that the speckles are multiply overlapped 
and both intensity and polarization fluctuations smear out. 
The polarization properties of radiation are determined 
by the second-rank polarization tensor. The correlation 
properties of the scintillation pattern are therefore de- 
scribed by a fourth-rank correlation matrix. If, however, 
the scattering medium is 'passive,' i.e., does not affect the 
polarization, then most of the information contained in 
this correlation matrix is redundant. The properties of ra- 
diation are then fully described by four scalar quantities, 
— the Stokes parameters, which are additive for indepen- 
dent (incoherent) radiation fluxes ( Rybicki fc LightmaE] 
19790 They are the intensity / and Q, C/, and V. The 
quantities Q and U describe linear polarization, and V de- 
scribes circular polarization. The degree of polarization is 
defined as 
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I^,,/I^{Q^ + U^ + VyiI 



(1) 



where /poi is the polarized intensity. Note that the polar- 
ized and total intensities are proportional to the integrated 
fluxes from the source which are directly measurable. 

We assume that scintillations do not affect polarization 
properties of radiation. This is valid if either the pro- 
jected turbulence is (on average) isotropic, i.e., there is 
no preferred direction in space, or the geometrical size 
of the inhomogeneities which contribute to scattering is 
much larger than a wavelength, s ^ A. Indeed, if s '--^ A 
and there is a preferred direction in space, then transmis- 
sion coefficients through a scintillation screen ('diffraction 
grid') of waves with different orientations of the plane of 
polarization with respect to that preferr ed direction will 



be different. It has been demonstrated ( Simonetti et al 
19840that birefringence in the interstellar medium is neg- 
ligible in discussion of scattering. Also, the effect of angu- 
lar Faraday depolarization across a sour ce is weak for the 
observing frequencies above ~ 10 MHz ( Linfield 1996 ). 

Let us assume that a point source is located at infin- 
ity. It illuminates a thin phase changing screen located at 
a distance D from an observer plane {x,y) and produces 
a speckle pattern on this plane. The statistics of speck- 
les is determined by the screen (i.e., the interstellar tur- 
bulence), and is usually characterized by the second mo- 
ment correlator of the complex electric field of the electro- 
magnetic wave, (i5(r)i?*(r + Ar)) = exp [— Z30(Ar)/2] ex 
exp [— (|Ar| /rdiff)^"^], where r = {x,y) and Ar are two- 
dimensional vectors, (...) = //... dxAy denotes a spatial 
average, the 'star' denotes complex conjugation, and (3 is 
the power-law index of the power spectrum of the electron 
density fluctuations, |(5ne((7)p oc q"^ with q being the spa- 
tial wave-number. The quantity D^ is the phase structure 
function which yields the root mean square phase shift 
along different paths. The value of /3 for the galactic in- 
terstellar medium is somewhat uncertain but close to the 



adopt the approximate value of /? « 4 for which the corre- 
lator is Gaussian, which greatly simplifies the calculation. 
The intensity distribution, Ip(x,y), in the speckle pat- 
tern produced in observer's plane by a point source is de- 
termined by the Fourier transform of the electric field cor- 
relator. We introduce the normalized intensity distribu- 
tion in the pattern and its Fourier transform as follows, 

W{x,y)^ip/{Ip) EE (/p- (/p))/(/p) ^ W(fc,,fc,), (2) 

where ^ denotes a Fourier conjugated pair and the 'calli- 
graphic' letters are used to denote the Fourier transformed 
quantities. Since the Stokes parameters are additive, the 
point source speckle patterns of Qp, Up, and Vp are also 
described by W. 

The intensity distribution produced by an extended 
source is obtained from the convolution of the intensity dis- 
tribution of a point source with the source surface bright- 
ness distribution, Pj(9x,9y), 

i{x,y) ^W[x,y)*Pi{9,,9y) 

= jj Wix -D9„y- D9y) Pi{9,,9y) d9., d9y 

^ yV{k^,ky) Vl{Dk^,Dky) = X{k„ky), (3) 

where {9x,9y) are angular coordinates in the sky, and 
we normalize the total flux to unity, i.e., (/) = 
JJ Pi{9x,9y)d9xd9y = 1. In equation (^), the convolu- 
tion theorem have been used. The amplitude of intensity 
fluctuations due to scintillations is determined by the mod- 
ulation index: 
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1/2 



(4) 



where the last expression follows from Parseval's theorem. 
By analogy with intensity, we define the distributions 
and modulation indices of the other Stokes parameters, 

Q^W^Pq, U:^W*Pu, V^W^Pv, (§) 

The amplitude of fluctuations of the polarized intensity is 
described by 



rriQuv = {mg + mfj + my) 



1/2 



(5) 



Kolmogorov theory prediction P — 11/3 (Armstrong et al 



By analogy with equation M), we formally introduce a 
measure of the 'degree' of polarization, m-^ = mquv /itii- 

3. ANALYTICAL MODEL 

An analytical calculation of the modulation indices is 
a difficult mathematical task, even for relatively simple 
models of source emission and interstellar turbulence. For 
this reason, we consider the simplest source model. We as- 
sume that the source has a uniform surface brightness over 
a rectangular angular region on the sky: —9s < 9^ < 9s, 
—9s < 9y < 9s, with 9s being the angular size ('radius') of 
the source. This approximation simplifies considerably the 
related integrals and provides a good approximation to a 
circular source. We normalize the total flux to unity. We 



choose the surface distribution of the Stokes parameters 
to be, 



1 
4l, 



■j^[ei0, + 9,)-Qi0,~0,)] 



X [Q{dy + 9,) ~ ei0y 
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^x, "y 



Pi 



U[0x,0y) 

Pv{9x,0y) 



' TTsPli 

0, 



r,0y] — cos — COS 



Tr0rr 



7r0„ 



sin — — sm — -p 
2 9s 2 9s 



(6a) 
,(6b) 

(6c) 
(6d) 



where tTs is the 'intrinsic' degree of polarization of radi- 
ation and Q{x) is the Heaviside step-function. The func- 
tions Pq and Pu are shown in Figure 0. These distribu- 
tions mimic the radiation field which islinearly polarized 
in the radial direction from the center of the source. We 
take the intensity distribution function to be Gaussian, 



Wix,y) 



V2i 



V^D0a 



exp 



DX 



(7) 



which is normalized so that {W"^)^/^ — nisc, where rrisc 
is the 'typical' modulation index for a point source. We 
do not specify tTs, ^O; and TOsc here and keep them as free 
parameters. The degree of polarization, tTs, depends on 
the emission mechanism and may be as high as ^ 70% for 
synchrotron radiation. The quantities ^q and nisc depend 
on the regime of scintillations and the wavelength of radi- 
ation. Namely, ^o ~ Tdis/D and TOsc — 1 for diffractive 
scintillations, and 0o ^ rrcf/D a nd ttIsc < 1 for refractive 
scintillations (see Narayan 1992 and also footnote 2). 

Using the formalism from the previous section (§0), 
we calculate the modulation indices mi, ruQuv, and ttt-jt. 
Upon lengthy and cumbersome mathematics, some details 
of which are presented in Appendix, we arrive at the fol- 
lowing expressions. 
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[B{d)-Cm' 



(8a) 
(8b) 
(8c) 



where ^ = ^/29s/0o is the size of the object in units of 
the characteristic speckle scale, and A{d), B{i}), and C{d) 
are complicated expressions involving the error functions 



of rea l and imaginary arguments, given by equations (A£)- 
( All ). Using equations (|4|),(^, it is now straightforward 
to calculate the quantities which are directly measured, 
namely the amplitude of fluctuations of total and polarized 
intensity (flux), m/ and mquv- The dependence of the 
scintillation indices on d is presented in Figure g. Clearly, 
for a small source size, (6*3 ^ 0q), there are no fluctua- 
tions of polarization, while the intensity fluctuations are 
maximum ~ 100% x Wgc- As the speckle size approaches 
the source size, the observed radiation is partially polar- 
ized. At the same time, the intensity contrast decreases 
due to the overlap between the speckles. The fluctuations 
of polarization peak when 0s '^ few x 9q at the value of 
~ 20% X TTs X nisc- For a larger source, the fluctuation 
amplitude of both polarized and total intensity decreases. 
However, the 'degree' of polarization (m^r = ^tiquv I'm^i) 
continuously increases with the increasing source size and 
asymptotes at ~ 70% x tTs in our model. Thus, the asymp- 
totic value of TO^ is nearly independent of the scatterer 
and yields the intrinsic degree of polarization at the source 
within a factor of order unity. 

4. CONCLUSION 

In this paper we propose a new observational tech- 
nique — the detection of polarization scintillations — to 
study the structure of the radiation field at compact ra- 
dio sources. We have developed a general formalism and 
demonstrated it on a simple analytical model. For ran- 
domly polarized radiation with the degree of polarization 
at the source typical of synchrotron radiation ~ 75%, the 
amplitude of fluctuations of the polarization flux is cal- 
culated to be ~ 15% in the case of diffractive interstellar 
scintillations and an order of magnitude lower in the case 
of refractive scintillations. 

We should note that these quantitative results were ob- 
tained for an over-simplified model and provide an illustra- 
tion of the predicted effect. Radiation fields of real sources 
are more complicated and are determined by the prop- 
erties of emission processes and radiative transfer. The 
statistical properties of the interstellar turbulence are de- 
termined by pumping and cascade mechanisms and the as- 
sumption that the fluctiations are gaussian is also rather 
crude. Therefore, we expect that numerical values of the 
modulation indices computed in a more detailed study may 
be different from those presented here. 

The author is grateful to Ramesh Narayan interest- 
ing and useful discussions and his comments on the 
manuscript and Avi Loeb for interesting discussions. This 
work was supported by NASA grant NAG 5-2837. 



APPENDIX 



CALCULATION OF THE SCINTILLATION INDICES 

First, we introduce new dimensionless variables, ^ = x/ D, rj — y / D. We use the following definition of the Fourier 
transform of W: 



Wik^,k^)^ JJ IU(C,77) e^p[~2m{k^^ + h,r,)] d^dr,, (Al) 

and similarly for other quantities. The Fourier transformed point source speckle distribution, equation (R), is 

W{k^, kri) = V2^m3c ^0 exp [-^^^^ ^^2 ^ ^2^j _ ^^2) 



The Fourier transformed distributions of intensity and the Stokes parameters, equations (p), expressed in terms of the 
new variables, are calculated to yield 

p,(fc^,fc^) = !i!i^!i!i^, (A3a) 



'PQih^^v) = T^s 



Trfcj 7rfc,j 

sin27rfc^0s /sin7r(2fc,,0s — 1) sin7r(2fc^6's + 1) \ sin27rfc,j6's /sin7r(2fcj6's — 1) sin 77(2/0^6*8 + 1 



2Tike 



7r(2fc,,6ls - 1) 7r(2fc,,6's + 1) 



27rfc„ 



7r(2fc^6ls - 1) 



Vu{k^,kn) = -TTs 



2 /sin7r(2A;56's - 1) sin7r(2fc^6's + 1)\ /sin7r(2fc^6's - 1) sin7r(2fc,,6's + 1) 



T:{2k^9s + 1) 

(A3b) 



T:{2k^es - 1) 



T:{2k^e, + 1) 



7r(2/c„6l, - 1) n{2k^es + 1) 



(A3c) 



The amplitude of fluctuations, which is the auto-correlation function at zero separation, may now be calculated using 
Parseval's theorem [see equation (Q)]: 



{P)^{i{C,v)ii^ + SLv + Sv)) 



55=0 

5ri = Q 



Wikc, kr,) W{-k^, -fc,,) Vl{k^,kr,) Vl{~k^, -kr,) dfcj dfc„ (A4) 



and similarly for (Q^) and (t/^). In calculating these quantities, one encounters the following basic integrals: 



h = 



_^2,2 sin TT^ 



^2 d$' -^2 = 



_^2^2 sin^ 7r(^ - 1) 



(^_1)2 d^' ^3= ' e 



_^2^2 sin^ 7r(C - 1) 



de 



(A5) 



Other integrals can be reduce to the above integrals by noticing that sin7r(^ — 1) = sin7r(^ + 1) = — sIutt^. To illustrate 
the method of taking these integrals, consider I2 as an example. We write I2 as a function of a new parameter, say a, so 
that it reduces to the original integral for a = 1, 



h(a) 



,^2^2 sin^ Q;7r(^ - 1) 



dC. 



Differentiating twice with respect to this new parameter, we obtain an integral which can now be evaluated, 



/2^^(a) = 27r2 / 6-""^' cos27ra(C - 1) d^ = 277^ ^ exp 
J-00 a 

One should now integrate l2aai'^) twice with respect to a and then set a = 1, 

^2 = / f / /2aa("i)dai j da. 



cos27ra. 



(A6) 



(A7) 



(A8) 



The lower limits in the integrals are taken so as to automatically include constants of integration by taking into account 
that /2(0) = ^2a(0) = 0, by definition. The integrals involving the error functions, which emerge after the first integration, 

may be taken by integrating by parts and using the fact that erf'(^) = (2/^/7F) e^'> . Finally, the integrals /1-/3 become, 
respectively, 

^3/2 



A{^) 



d 



1 + V7ri9erf(z9) 



B{d) = 



r3/2 



1 „ V^p-^V'S' 

2 



G-^)"«e-^)+G+'^)"«(^'^)-r*© 
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erfi( '--id]+ crfi( ^ + i ^"j - 2 erfif^ 

2 ,2 



(A9) 
(AlO) 
(All) 



where erf(x) = 2/ ^ j^ e * dt is the error function, erfi(z) — 2/y^/p e* At — — ierf(i z) is the error function of imaginary 
argument, and d — \/2^s/^o- The scintillation indices are then readily calculated to yield equations (ph. 
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Fig. 1. — Distributions of the Stokes parameters, Pq and Pjj, over the source in the analytical model. From black to white, P varies from 
-1 to +1. 

Fig. 2. — Modulation indices of the intensity (m,/), polarized intensity (mQijv), and 'degree' of polarization {m,r = 'tiq(7v/"1/), as 
functions of the normalized source size i? for completely polarized radiation, tts = 1, in the diffractive scattering regime, rrisc = 1- The thin 
vertical line marks where 6= = 6n. 
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